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TECHNICAL MEMORANDUM X-64984 


FEASIBILITY DEMONSTRATION 'FOR HYDROGEN CHLORIDE 
DETECTION USING A CHEMISORPTION TECHNIQUE AND 
A QUARTZ CRYSTAL MICROBALANCE 


I. INTRODUCTION 


There is a growing awareness and concern about the atmospheric pollu- 
tant gases tli at are being distributed in our environment. One of the gases of 
concern is hydrogen chloride (HCl) , which is extremely hydroscopic. The 
concerns range from the effects gaseous and/or aqueous concentrations may 
have on individual health to the effects of reactions that contribute to the deple- 
tion of tlie ozone layer. 

A letter from the National Research Council expresses the limits to 
which an individual may be exposed as 8 parts per million (ppm) and further 
states: 


"The STPL’s (Short-term Public Limits) for HCl were established 
by the committee on Toxicology in August 1971. They are: 

10 minutes at 4 ppm 
30 minutes at 2 ppm 
GO minutes at 2 ppm 
1 hour daily at 2 ppm 

5 hours/day, 3 to 4 days /month at 0.7 ppm 

• • • 

The other PEL’S (Public Emergency Limits) for HCl are: 

30 minutes at 20 ppm 
60 minutes at 10 ppm." 

Therefore, it is necessary that suspected high concentrations or 
repeated exposures to HCl in the troposphere be measured to insure personal 
safety. 

Another concern is expressed by Dr. Fred Kaufman, a chemical kineticist 
at the University of Pittsburgh, and his research group who are studying the 
reactions involving HCl which lead to the depletion of the ozone layer, such as 
the sequence, 

HCl + OH — » H 2 0 + Cl 

Cl + O 3 — CIO + o 2 . 


These reactions result in a net loss of ozone molecules which may con- 
tribute to depletion of the ozone layer. 

While many methods and instruments are available to measure the con- 
centration of HC1, most of them are expensive, cumbersome, and delicate. An 
inexpensive, portable, and rugged device which can be refurbished easily and 
remotely operated would be extremely valuable. Such a device has been con- 
ceived and its feasibility demonstrated. The device combines chemisorption 
techniques with a quartz crystal microbalance. The theory, application, and 
demonstration of this device are discussed in the remainder of this report. The 
device is called a CATCH (continuous accumulation and time correlated history) 
detector. 


II. THEORY 


Basically, we must establish a theoretical model and relationship between 
the measurement parameters and the gas we desire to monitor. To accomplish 
this we will initially summarize the Maxwell-Boltzmann law in classical statis- 
tical thermodynamics, adapting it to die microbalance measurements of HC1. 

It mil be noted that chemisorption is assumed here to be time independent. 

The molecular model employed is based on the following six classical 
assumptions: 1 * 


1. ANY MACROSCOPIC VOLUME OF A GAS CONTAINS A 
VERY LARGE NUMBER OF MOLECULES. This assump- 
tion is justified by all experimental evidence. The number 
of molecules in akilomole (Avogardro’s number N ) is 

A 

6. 03 x 10 2G . . . . At standard conditions, 1 kilomole of an 
ideal gas occupies 22. 4 m 3 . Hence at standard conditions 
there are approximately 3 x 10 20 molecules in a cubic 
meter, 3 x 10 19 in a cubic centimeter, and 3 x 10 1G in a 
cubic millimeter. 

2. THE MOLECULES ARE SEPARATED BY DISTANCES THAT 
ARE LARGE COMPARED WITH THEIR OWN DIMENSIONS 
AND ARE IN A STATE OF CONTINUOUS MOTION. The 


1. Francis W. Sears and Gerhard L. Salinger: Thermodynamics, Kinetic Theory, 

and Statistical Thermodynamics, Third Edition, Copyright 1975, Addison- 

Wesley Publishing Company, Reading, Massachusetts, pp. 251-252. 
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diameter of a molecule, considered as a sphere, is about 2 
or 3 x 10“ 10 m. If we imagine one molal volume at standard 
conditions to be divided into cubical cells with one molecule 
per cell, the volume of each cell is 30 x 10 -27 m 3 . The 
length of one side of such a cell is about 3 x 10” 9 m, which 
means hi at the distance between molecules is of this order 
of magnitude, about 10 times the molecular diameter. 

3. TO A FIRST APPROXIMATION, WE ASSUME THAT 
MOLECULES EXERT NO FORCES ON ONE ANOTHER 
EXCE PT WH EN THEY COLLIDE . Therefore between 
collisions with other molecules or with the walls of (lie 
container, and in the absence of external forces, they 
move in straight lines. 

4. COLLISIONS OF MOLECULES WITH ONE ANOTHER AND 
WITH THE WALLS ARE PERFECTLY ELASTIC. The 
walls of a container can be considered perfectly smooth, 
so that there is no change in tangential velocity in a 
collision with the walls. 

5. IN THE ABSENCE OF EXTERNAL FORCES, THE MOLE- 
CULES ARE DISTRIBUTED UNIFORMLY THROUGHOUT 
THE CONTAINER. If N represents the total number of 
molecules in a container of volume V, the average number 
of molecules per unit volume, n, is 

n = N/V . 


The assumption of uniform distribution then implies that in 
any element of volume A V, wherever located, the number 
of molecules AN is 


AN = nAV . 


Obviously, the equation above is not correct if AV is too 
small, since the number of molecules N, although large, 
is finite, and one can certainly imagine a volume element 


so small that it contains no molecules, in contradiction to 
the equation above. However, it is possible to divide a 
container into volume elements large enough so that the 
number of molecules per unit volume within them does not 
differ appreciably from the average, and at the same time 
small enough compared with the dimensions of physical 
apparatus that they can be treated as infinitesimal in the 
mathematical sense and the methods of differential and 
integral calculus can be applied to them. For example, a 
cube l/lOOO mm on a side is certainly small in comparison 
with, the volume of most laboratory apparatus, yet at 
standard conditions it contains approximately 30 million 
molecules. 

6. THE DIRECTIONS OF MOLECULAR VELOCITIES ARE 
ASSUMED TO BE DISTRIBUTED UNIFORMLY. To put 
this assumption in analytic form, imagine that there is 
attached to each molecule a vector representing the 
magnitude and direction of its velocity. Let us transfer 
all these vectors to a common origin and construct a 
sphere or arbitrary radius r with center at the origin. 

The velocity vectors, prolonged if necessary, intersect 
the surface of the sphere in as many points as there are 
molecules and the assumption of uniform distribution in 
direction means that these points are distributed uni- 
formly over the surface of the sphere. The average 
number of these points per unit is 

N 

47a’ 5 * 


and the number in any element of area AA is 

AN =i^ AA ■ 


wherever the element is located. As in the preceding- 
paragraph, the area must be large enough (that is, it 
must include a large enough range of directions) so 


that the surface density of points within it does not differ 
appreciably from the average. Because of die large 
number of molecules, die range of directions can be made 
very small and still include a large number of points. 

Using polar cooi’dinates (r, 0, <p ) , the number density of molecules witii 
velocities having directions between 0 and 0 + A0 , and 0 and 0 + A 0 is 


An = -j- sin 0 A 0 A 0 . (l) 

00 4ir w 

Not all molecules will have die same speed, but ratiier will have speeds tiiat 
range from zero to infinity. 

It is assumed that in an equilibrium condition die number of molecules 
witii velocity in any specified range will remain constant; i. e. , altiiough a given 
molecule may lose or gain speed and change direction due to a collision with 
another molecule and/or a surface, some other molecule in the system is also 
changing its velocity so that the distribution of velocities remains constant. 

The molecular flux (<T>) at any surface is defined as the total number of 
molecules arriving (AN) at the defined area (A A) ner unit time (At) ; hence, 


AN 

AA At 


( 2 ) 


The construction of an arbitrary plane surface in space and the development of 
an equation for the molecules that will be incident on that surface results 
in the classical formula for the flux: 


4> - 1/4 n v 


( 3 ) 


where n is the number of molecules/m 3 , and v is die average velocity of 
die molecules in m/s. The flux $ is tiien the number of molecules/m 2 /s. 


Using the Maxwell- Boltzman speed distribution function 



the definition of the average or arithmetic mean speed 

v = i 


and approximating the sum by an integral we get 



The ideal gas law gives us 

PV = NkT or f = ur ^ n ( 6 ) 

Then by substituting equation ( 5) into equation (3) , we obtain the Maxwell- 
Boltzmann flux which can be expressed in terms of the number of molecules/m 3 
(n), the temperature (T) in K, the mass (m) in kg, and Boltzmann’s constant 
(k) = 1. 3806 X 10~ 23 J/K; 



In equation (6) , the number of molecules/m 3 can be expressed in terms 
of the pressure (P) in N/m 2 , Boltzmann's constant (k) , and the temperature 
(T) in K. Theoretically, it can be calculated from equation (6) that the number 
of molecules/m 3 at standard temperature (273 K) and pressure (1. 0132 X 10 5 
N/m 2 ) is; 



2. 7 x 10 25 molecules/m 3 or 2. 7 x 10 19 molecules/cm 3 . 
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For air quality measurements, we need to relate parts per million or 
billion, the accepted terms for air quality discussion, to molecules/cubic 
centimeter. Therefore, if we assume that there is 1 ppb (l ppb = lx 10" 9 ) of 
atmospheric pollutant gas at standard temperature and pressure (STP) , then 
the number of these molecules (n*) will Joe one billionth the total number of 
molecules, or 


n’ = " = 2. 7 x 10 10 molecules/cm 3 


To determine the flux, we must identify the atmospheric pollutant gas 
so that we may obtain either the average velocity or mass to use in flux equation 
(7). For the purpose of this report, we will identify HCl. Since HC1 is hydro- 
scopic, the probability that the HCl molecule will attach itself to a water mole- 
cule must be taken into account. For the sake of simplicity and to be on the 
conservative side, let us assume that each and every HCl molecule has an H ? D 
molecule attached to it. The mass (m f ) of such a combined molecule is 9. 04 x 
10" 26 kg/molecule, and the average velocity, at a temperature of 273°K, calcu- 
lated from equation (5) would be 325. 8 m/s or 325. 8 x 10 2 cm/s. 

Substituting these results into equation (7) to obtain the flux, when the 
concentration level is 1 ppb, yields 

4> = -|n f v = 2. 2 x 10 14 molecules/cm 2 /s 

Assuming that a surface with a chemisorption coefficient of unity is 
available for these incident, combination molecules, the rate of accumulation 
would be 


4> m* = 1. 99 X 10 ~ 8 gm/cm 2 /s 

This is substantiated by an alternative approach using de Boer’s equa- 
tions which give 

$ m’ = (3. 52 x 10 22 ) = 1. 9 X 10“ 8 gm/cm 2 /s 

(mT) /z 

where p- (l ppb) 760 x 10“ 9 mm of mercury, T = 298°K, and m = 54. 5 gras. 
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Therefore, it is possible to measure low concentrations of atmospheric pollutant 
gases provided die relationship between the sensitivity of the detector and the 
adherence of the molecules to the surface of die detector is correctly balanced; 
i. e. , if die chemisorption coefficient, defined as die ratio of molecules of 
interest that adhere to the surface divided by die total number tiiat come in 
contact witii it, is 1 and die sensitivity of die detector is 2 x 1C" 8 gm/em 2 /s, 
tiien die concentration of 1 ppb can be measured as discussed. If the chemisorp- 
tion coefficient is 0. 1 and die sensitivity rf die detector is 2 x 10" 9 gm/cm 2 /s, 
die concentration of 1 ppb would again be measurable. Therefore, it is die 
relationship between die chemisorption coefficient and the sensitivity of die 
detector tiiat determines die measui’able concentration. 


til. APPLICATION 


The detection of an atmospheric pollutant gas like HC1 in low concentra- 
tions car. be accomplished as discussed in die preceding section. The problem 
is broken into two parts: die sensitivity of die detector and die adherence of die 
molecules to die detector surface. For tiiis discussion, assume die chemisorp- 
tion process to be independent of time. 

Since there exist a multitude of techniques and microbalances that cover 
the sensitivity scale do mi to 10" 9 gm/cm 2 , this aspect of the problem merely 
becomes one of selecting the proper means for measuring die range desired. 

The real-life problems make it a little more difficult than implied in die above 
statement; however, in practice it is still a process of choice. 

The second aspect, developing or devising a method whereby die mole- 
cules of interest can be coupled to die measurement surface, remains as a 
problem to be solved. 

For the purpose of this report and for a more detailed discussion, we 
have identified HC1 as a gas of interest and have assumed a concentration level 
of 1 ppb. Because this is a very small molecule, even with die water molecule 
attached, it is very lightweight, and a very sensitive microbalance is required. 

In tiiis case, let us choose a 10 MHz quartz crystal microbalance whose 
sensitivity is on the order of 10" 9 gm/cm 2 -Hz. The sensitivity of a particular 
quartz crystal is found by using die formula 2 


2. A. W. Warner: Micro Weighing witii Quartz Crystal Oscillators — Theory and 
Design, Chapter 5, p. 143. 
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m”/(AAf) = (1/f) p t 


( 8 ) 


where m” is the mass collected in grams, A is the collecting area in cm 2 , 

Af is the change in frequency (Hz) , f is the crystal frequency (Hz) , p is the 
density of the quartz in gm/cm 3 , and t is the thickness of the crystal in cm, 
Assuming f- lx 10 7 Hz, p- 2.56 gm/cm 3 , andt= 0.016 cm, the sensitivity 
(s) would be 4. 2 x 10" 9 gm/cm 2 /Hz. 

To simplify calculations, we will derive a formula for the concentration 
level (c) as a function of the parameters of the molecules of interest and the 
quartz crystal microbalance: 


number of molecules of interest 
total number of all molecules 


n’ 4<&kT 



where 


= (molecules/m z /s) = 


( 9 ) 



and 


s' = (lx lO" 3 kg/gm)(lx 10 4 cm 2 /m 2 )s = 4. 2 x 10 “ 8 kg/m 2 /Hz ; 
therefore, substituting equation ( 10) into equation ( 9) yields 


4( Af/At) s’lcT / A f s ’ \ / 1 \ 4kT 

/ 8 kT\V 2 \ At A m V P 


(ii) 


where Af is the change in frequency (Hz) , At is the elapsed time (s) , s’ is 
the sensitivity of the quartz crystal (kg/m 2 /s) , k is Boltzmann’s constant 
(J/K), T is the temperature (K) , m is the mass of the molecule (kg) , P is 
Hie pressure (N/m z ), and v is the average velocity of the molecules (m/s) , 

With this formula, we can solve for the concentration, assuming that all 
the molecules of interest, and only these molecules, are adhering to the surface 
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of tile quartz crystal. To have confidence in tills assumption it is necessary to 
provide a surface which is selective and which has a known stay time or sticking 
coefficient. 

Athens College, working through contract NAS8-28058, recommended 
that an aromatic amine [triphenyl amine f EC c H 6 J 3 ) j be used because tertiary 
amines are less reactive toward water vapor and H 2 S0 4 if the substituents ( in 
this case, 3 phenyl groups) are large and nonpolar. The low vapor pressure 
of large molecular weight amines also increases their usefulness for mass 
measurements with the quartz crystal microbalance. 

The successful technique for bonding this amine to the surface of the 
quartz crystal was to vapor deposit the triphenylamine on the crystal surface in 
a vacuum chamber. The amount of triphenylamine can be monitored during 
this process, and reproducible coatings are possible. Since the boiling point of 
triphenylamine is approximately 365°C, the coated CATCH detector is 
durable to the high temperature limit of the quartz crystal microbalance. 

Several quartz crystals of the sensitivity described earlier (4. 2x 10 -9 
gm/cmVs) were coated with the amine and exposed to HC1 and other mixtures 
and substances. While these preliminary tests were not calibrated, they did 
show the selective nature and some of the characteristics between the amine, 
HC1, and the quartz crystal microbalance. 

At this point, a chance to demonstrate the feasibility of such a device in 
the field by measuring HC1 as a by-product of solid rocket propellant combustion 
during Viking launches from Kennedy Space Center was provided through the 
efforts of Dr. Briscoe Stephens. He originally stimulated the investigation into 
tiie possible application of such a device in the measurement of HC1 while he 
was working with the Joint MSFC-LaRC-KSC (Marshall Space Flight Center- 
Langley Research Center-Kennedy Space Center) Rocket Exhaust Effluent 
Diffusion prediction and Monitoring Team. Because the Viking launches from 
KSC represented the last potential opportunity to field test this "new” technique 
in conjunction with other supporting instruments under actual launch conditions, 
further calibration was postponed. However, it is recognized that the chemi- 
sorption coefficients and parameters are needed together with a more detailed 
analysis of the pollutant gas reaction characteristics. 


IV. DEMONSTRATION 


The opportunity to demonstrate feasibility of flic CATCH detector was 
provided at the Viking B launch from Cape Kennedy on August 20, 1975. It 
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occurred before die calibration of die instrument could be completed. However, 
since it represented die last opportunity to field test this new technique, five 
crystals were coated and dispersed at three locations around the launch site. 
Only one of these three locations, with two operating detectors, was on a line 
from die launch site over which a portion of the exhaust cloud passes at a height 
of several hundred meters. This location was 4 km from the launch site, and 
die data recorded are presented in Figure 1. 

Several facts are necessary to interpret the results: 

1. Data points arc 3 minutes apart due to the time-sharing sequence of 
the recording instrumentation. 

2. Data points are the average of 44 frequency readings taken during a 
1-minute interval. The standard deviation (a) is 50 Hz. 

3. An increasing frequency indicates a decrease in the temperature. 

4. A decreasing frequency indicates an increase in the thickness or 
mass addition to die crystal. Therefore, detection of HC1 is a negative-going 
peak. 


5. The two detectors were at different heights above ground level. 
Detector A was 4 m (13 ft) above the ground and detector C was 1. 8 m (0 ft) 
above the ground. 

Keeping these facts in mind, the data in Figure 1 indicate an abrupt 
change in the temperature of detector C about 5:35 p.m. A similar, but smaller, 
shift occurs in detector A. Unfortunately, die temperature at this location was 
not monitored; therefore, a quantitative verification cannot be made. However, 
a severe thunderstorm moved in from the southwest and blocked out the sun, 
causing an abrupt, physically noticeable change in the temperature as noted by 
the operator during this time period. During the other time periods the rate of 
increase in frequency is regular, which would be indicative of a slowly changing 
temperature. 

Between 5: 40 and 5: 50 p. rn. both detectors record a negative-going 
peak of the same amplitude ( ~ 150 Hz) . Both, peaks occur during the time the 
exhaust cloud was seen visibly to pass overhead as recorded in the log book by 
the operator. The lowest portion of the visible exhaust cloud was several 
hundred meters above the ground. The amplitude of the peaks is 3 a and, there- 
fore, from statistics represents less than a 5 percent; probability of being a 
random deviation. For these reasons, it is concluded that these events repre- 
sent a measurement of HC1 with a rise time of 180 s. Using these values in 


11 


FREQUENCY MINUS A BIAS (HERTZ) 







equation <11) for the concentration level, the calculated concentration <C c ) at 
STiP is 0. 2 ppb. 

Preliminary laboratory experiments indicate that the ratio of incident 
molecules to the number that are chemisorbed is not unity. Assuming a low 
.value of 0.001, as suggested by the laboratory experiments, for the chemisorption 
ratio (a) gives 

C = — C = 200 ppb or 0. 2 ppm . 
tv c 


The data also show smaller negative peaks later which are not designated 
as HC1 measurements. While it is possible that they are HCl indicators, either 
they are not detected by both CATCH detectors at the same time or their ampli- 
tude is less than 2a. Therefore, they are left unexplained. 

It can be seen from Figure 1 that the capture and measurement process 
was followed by a loss of the HCl molecules. This phenomenon is in contra- 
diction to the premise of chemisorption upon which the theory is based. The 
laboratory tests indicate that this effect is real and consistent with measure- 
ments of HC1 in the atmosphere using this technique. Investigations are 
presently being pursued to establish the basis for the phenomenon and why it 
does not follow the theorized chemisorption process. 


V. CONCLUSION 


There is a growing awareness and concern for atmospheric pollutant 
gases which are being distributed in our environment. IIC1, one of these gases, 
causes concerns which range f rom the effect gaseous and/or aqueous concen- 
trations rnay have on individual health to the effect of reactions that contribute 
to the depletion of the ozone layer. 

It is, therefore, important that an HCl monitoring system be pi’oposed 
and demonstrated that would be small, low cost, portable, and remotely operable. 
These features would provide a means of measuring concentrations of HCl in 
locations that vary from the open troposphere to the freshman chemistry 
laboratory. 

A device which meets these criteria has been described and discussed 
in this document. It couples a chemisorption technique with a real-time, mass- 
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monitoring quartz crystal microbalance. The quartz crystal microbalance is 
sensitive to approximately 10" 9 gm (10 13 molecules of HC1) and is, therefore, 
capable of measuring concentrations of HC1 down to 1 ppb at STP. The feasi- 
bility of this device was demonstrated at the Viking B launch at Cape Kennedy 
on August 20, 1975, where HC1 was a by-product of the solid propellant com- 
bustion. The data from this test indicated a concentration level of 0. 2 ppm at 
4 km from the launch site based on a chemisorption coefficient of 0.001. 
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